Nuclear spin polarization in a single quantum dot pumped by two laser beams 
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We theoretically investigate dynamic nuclear spin polarization in a self-assembled quantum dot 
pumped optically by two laser beams. With the assumption that a noncollinear interaction between 
the hole spin and nuclear spins leads to nuclear spin polarization, we find that both weak and 
strong nuclear spin polarizations can arise, depending on the intensities and central frequencies of 
the lasers. For weak nuclear spin polarization, we use a perturbation method to show that the 
distribution of the nuclear spin Overhauser field may become significantly narrower. Using Monte 
Carlo simulations to study a single quantum dot, we find that strong nuclear spin polarization can 
also be generated via appropriate optical pumping. 

PACS numbers: 73.21. La, 78.67.Hc, 72.25.Fc 



I. INTRODUCTION 

Dynamic nuclear spin polarization by optical pump- 
ing on single quantum dots (QDs) is extensively studied 
in recent years^— due to its importance for improving 
techniques in quantum computing, and for understanding 
the nuclear spin environment of the electron (hole) spin 
qubit. For instance, an efficient suppression of nuclear 
spin fluctuation by optical pumping allows one to pro- 
long the electron spin coherence time in a QD i 13 ' 25 Also, 
significant nuclear spin polarization may be achieve d 2 ' 3 ' 11 
even by optical excitation of spin-forbidden transitions^ 
In order to understand the dynamic nuclear spin polar- 
ization induced by continuous optical pumping on single 
QDs, a microscopic theor y 31 ' 32 was recently introduced 
to explain relevant experiments! 1 ' 26 However, this the- 
ory is restricted to single laser pumping and only treats 
weak nuclear spin polarization, while experiments involv- 
ing two laser beam o 13 ' 25 or strong nuclear spin polariza- 
tion remain to be explained microscopically. 

In this paper, we study dynamic nuclear spin po- 
larization in a single QD charged with one electron 
and pumped by two narrow-linewidth continuous wave 
lasers (see Fig. [1| . We find that both weak and strong 
nuclear spin polarizations can be generated. Specifically, 
when the laser with central frequency u)\(p\ and moder- 
ate Rabi frequency f^i(2) is off-resonant (resonant) with 
the transition between the exciton and the electron spin 
eigenstate \x — (£+)), the nuclear spins tend to have 
weak polarization for ft i <C 0,2 , but strong polariza- 
tion for fix 3> ■ In the former case, we derive a 
Fokker-Planck equation^ 3 - for the evolution of the prob- 
ability density of nuclear spin polarization. We use the 
Fokker-Planck equation to show that the nuclear spin 
fluctuation can be reduced, thus enhancing the electron 
spin coherence. Also, we perform numerical study of 
several interesting phenomena observed in laser spec- 
troscopy experiments f^ 3 - For the case of strong nuclear 
spin polarization, the perturbation method for deriving 
the Fokker-Planck equation breaks down. We then cre- 
ate a toy model comprising a small QD in order to obtain 
unbiased results using the Green's function Monte Carlo 



simulation! 34 ' 35 Our Monte Carlo simulation shows that 
at least as high as 50% of the nuclear spin polarization 
degree can be generated by pumping optically on the sin- 
gle QD in our toy model. We note that a 50% degree of 
nuclear spin polarization is close to the large nuclear spin 
polarization degree observed in experiments by optically 
pumping single QDs^ i 11 ' 20 ' 36 

The paper is organized as follows. Section II briefly in- 
troduces the system Hamiltonian. Section III solves the 
equation of motion for the nuclear spin population. In 
Sec. IV, we study the case corresponding to the experi- 
ments in Ref. Il3l . where only weak nuclear spin polariza- 
tion was observed. In Sec. V, we first identify the factors 
responsible for large nuclear spin polarization. Then we 
numerically study a small QD of artificial size to show 
that large nuclear spin polarization can indeed be gener- 
ated by optical pumping on single QDs. We conclude in 
Sec. VI. 



II. MODEL 

Consider a self-assembled InAs QD charged with one 
electron, where the confined electron interacts with N 
indium nuclear spins and N arsenic nuclear spins. Defin- 
ing z as the growth direction of the QD, we apply an 
in-plane static magnetic field B in the x direction and 
label the electron spin eigenstates as |x±), as shown in 
Fig. [TJ Two linearly polarized coherent laser beams, one 
with the polarization of the electric vector along x di- 
rection, and the other along the y direction, selectively 
couple the two electron spin states to a common trion 
level, denoted as |T— ). The trion consists one heavy 
hole and two electrons^ with the two electrons in the 
singlet configuration. The Hamiltonian of the system is 



H = H r 



(1) 



where H c h is the Hamiltonian of the electron-hole system, 
H n the nuclear spin Zccman term, and H H1 the hyperfine 
interaction between the electron/hole spin and the nu- 
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where a (h).j is the hyperfine interaction strength be- 
tween the electron (hole) spin S e (h) and the nuclear spin 
Ij, the superscript x, y or z denotes component of the 
spin moment in the corresponding direction, and f3 = 
\P\e lS is the heavy- light hole mixing coefficient <22r— Be- 
sides the dipole-dipole coupling in the hole spin-nuclear 
spin interaction, there are extra hyperfine couplings with 
strength proportional to \fi\ in Eq. ((SJ that arise from 
the mixing between the heavy and light hole bands. The 
band-mixing is caused by an in-plane strain of the QD 
when one grows QDs in the Stranski-Krastanow growth 
mode during molecular beam epitaxy 40 ' 43 Observed val- 
ues of |/3| for typical In(Ga)As QDs range from 0.02 to 
0.7i 13 i 39 ' 44 i 45 The angle 5 is determined by the strain 
detai l 40 ' 43 of the QD and can in principle take any value 
between and 2ir. For simplicity, we have ignored the in- 
trinsic interactions (including dipole-dipole interactions) 
between nuclear spins in H n . The latter will be consid- 
ered later when studying relevant experiments, involving 
nuclear spin depolarization^ 

For study involving small nuclear spin polariza- 
tions, we consider a QD containing N — 9500 InAs 
molecules (estimated from Ref. [l3|). Without loss of gen- 
erality, we assume g e ^ = 0.49 (—0.13). Although this 
choice of parameters is made in order to study the experi- 
ments in Ref. [13, we note that our general results for both 
small and large nuclear spin polarizations do not rely on 
these specific assumptions. In the electron-hole system 
shown schematically in Fig. [TJ the laser connecting \x— ) 
and \T— ) may be slightly off- resonant with the transi- 
tion between the electronic and excitonic states, while 
the laser connecting |x+) and \T— ) is resonant unless 
otherwise specified. For this case with a strong external 
magnetic field of \B\ = 2.64 T, we can prove that the 
level \T+) decouples from the rest of the system (Ap- 
pendix |XJ, resulting in a three-level electron- hole system 
which we call a A system (AS). 



III. NUCLEAR SPIN DYNAMICS 

Following Refs. [H, and H3, we assume that the 
effective noncollinear hyperfine interactions between the 
electron/hole spin and nuclear spins can lead to dynamic 
nuclear spin polarization in a single QD. For simplicity, 
we only consider the hole spin-nuclear spin noncollinear 
interaction of the form S^Ij here, which appeared in 
Eq. @. We neglect the influence of the transverse part of 
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FIG. 1. (Color online) Schematic of an electron-hole system. 
The x-polarized laser connecting \x—) and \T— ) has central 
frequency cji and Rabi frequency fii , and the j/-polarized laser 
connecting \x+) and \T—) has central frequency u>2 and Rabi 
frequency Q, 2 . 



the electron (hole) spin-nuclear spin hyperfine interaction 
on the nuclear spin dynamics as discussed in Appendix [A1 
Using the theory developed in Ref. l3ll we derive an 
equation for the time evolution of the nuclear spin popu- 
lation P(t), i.e., the diagonal part of the reduced density 
matrix of the nuclear spins, 



Jt P(t)*--£{[ir, itw^P(t) 



[V 



l7W a ..-Pl 



(*)]}, 



(3) 



where ay denotes In (As) when there is an indium (ar- 
senic) atom at site j, and 
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where /x N is the nuclear magneton, p^ ss ' is the reduced 
density matrix for the steady state of the AS,— Tr de- 
notes trace over the A-systcm degrees of freedom, T is 
the time ordering operator, iJ e h,r(i') is obtained by per- 
forming a rotating frame transformation on iJ c h (see Ap- 
pendix [AJ, and the integration is from —t to t. However, 
since the scale of t is much larger than the time scale of 
the AS, we can perform the integration from — oo to oo 
and the final integration can be calculated through the 
quantum regression theorem] 48 ' 49 

There is a feedback loop between the nuclear spins 
and the AS. First, the rate W aj ,± in Eq. © for flip- 
ping the nuclear spins depends on the steady state of 
the AS. Second, the steady state of the AS depends on 
the nuclear spin state, since the eigenenergics of the two 
electron spin eigenstatcs |x±) in the AS are shifted by 
an effective magnetic field contributed by a e-jlj m 

H HI . Here we neglect the nuclear spin Overhauser field 
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contributed from the hole spin-nuclear spin hyperfine in- 
teraction since |a n .,-| is much smalle r 21 ' 43 ' 50 ' 51 than \a e j\. 
The effective magnetic field, i.e., the nuclear spin Over- 
hauser field, can be written as [see Eq. (|A7I) ] 



ft = [-Ai n is In + Aa s Jsas], 



(5) 



where 7 = 9/2, J = 3/2, *4i n (As) is the hyperfine constant 
of the indium (arsenic) nuclear spins^ and s In , s As are 
the expectation values of 
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Because of the nuclear spin Overhauser field ft, the nom- 
inal detunings 

Ai,o = Eif — — u>i, A2,o = Ej! — E x — U2 (7) 



are replaced by the actual detunings 



Ai = A li0 
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in the optical Bloch equation [see Eq. (|A15[) ]. Here 
E x , Ex and E^ are the respective eigenenergies (when 
there is no net nuclear spin Overhauser field) of \x), \x—) 
and |T— ), and Wx(2) is the central frequency of the laser 
beam that couples |x=p) and |T— ). For different adiabatic 
nuclear spin Overhauser fields on the AS, the Zeeman en- 
ergy of the electron spin is effectively shifted by different 
amounts, giving different flipping rates in Eq. ([3]). When 
h keeps on changing, the flipping rates in Eq. (|3]) keep 
on changing. This feedback loop between the AS and 
the nuclear spins in the QD will not cease whenever the 
distribution of ft is changing. 



IV. SMALL NUCLEAR SPIN POLARIZATION 



A. Approximation when |s In 



(As) 



< 1 



In practice, it is not easy to solve the distribution of ft 
from Eq. since the number of nuclear spins is large. 
Nonetheless, when the polarizations of the two species 
of nuclear spins are much smaller than 1, we can obtain 
solution for the evolution of the joint probability density 
of the nuclear spin polarizations 



p s (si„,s As ) = Tr P(t)5 hn!Sln 5g As , SAs , 



(9) 
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in(As) i s the delta function. In particular, 



we derive a Fokker-Planck equation 
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where v a and D a , given by Eq. (|P311|) . are known as the 
drift and diffusion coefficients^ 3 - respectively. From the 
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FIG. 2. (Color online) Flipping rate r In (A s ), P 



as a func- 



tion of the nuclear spin Overhauser field ft in a QD with 
N = 9500 InAs molecules. Parameters: B = -2.64 T, 
g e(h) = 0.49(— 0.13), |/3| = 0.2,5 = tt/9, a CJ = -10a htj = 
Aj/N, and 7s = 0.06, Ai, = A 2 ,o = 0,T = 0.4, fi t = 
0.24,^2 = 1-35 (unit: GHz). Here the gyromagnetic ratio 
is gjfin = 0.0093(0.0073) GHz/T, and the hyperfine constant 
Aj — 13.6 (U.l) GHz (taken from Ref . l52h when there is an 
indium (arsenic) nuclear spin at j. 



relation between s In , s As and h in Eq. ([5]), we can obtain 
the distribution of h from p s (s ln , s Aa ). 

In order to study the absorption coefficient of the probe 
laser in laser spectroscopy experiments on single QDsji 3 - 
we consider the dynamics of the mean of the nuclear spin 
Overhauser field. To simplify the calculation, we neglect 
the finite width of the distribution of the polarization 
s In and s As . Then the mean of nuclear spin polarization 
Sin(A S ) obeys (Appendix [B| 
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and similarly for v As . Here, (• • • ) denotes the respec- 
tive expectation value. In Eq. (|11[) we have added a nu- 
clear spin depolarization channel with rate 7d ep _ 0. The 
reason for adding this term is that, in the nuclear spin 
Hamiltonian H n , we have ignored the direct dipole-dipolc 
interactions which play an important role in depolarizing 
the nuclear spins4£ A similar depolarization term was 
also used in the numerical fitting of the experimental re- 
sults in Ref. [l3|. By choosing the electron spin dephasing 
rate 7 S = 0.06 GHz (this value is taken from the sup- 
plementary material of Ref. [l3l ). the rate of spontaneous 
decay from \T— ) to either \x+) or \x—) as T = 0.4 GHz, 
the Rabi frequency as f?i(2) = 0.24 (1.35) GHz for the 
probe (pump) laser beam, we calculate the flipping rate 
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FIG. 3. (Color online) Solid (dashed) curves show the steady 
distributions of the nuclear spin Overhauser field h built up 
with (without) optical pumping. Ai,o = —0.3,0,0.675 and 
±2 GHz in (a), (b), (c) and (d), respectively. Other parame- 
ters are the same with Fig. [2] except Ai n = 14.2 GHz here. 



r Tn (As),p as a function of h for A^o = A2.0 = 0, and 
obtain the results shown in Fig. [2] 



B. Numerical results 

There are two important observations made in Ref . Il3l . 
The first is an enhanced electron spin coherence time, 
which was observed indirectly via a deeper dark state 
dip of the probe laser absorption. The second observa- 
tion concerns a series of phenomena caused by dynamic 
nuclear spin polarization, which were observed via the 
change of the probe laser absorption in reference to the 
expected absorption when there is no dynamic nuclear 
spin polarization. We will study these two types of phe- 
nomena. 

For the first observation where optical pumping re- 
duces the nuclear spin fluctuation, we will study the 
steady distribution of the nuclear spin Overhauser field. 

For the second observation, we will simulate the time 
evolution of the absorption coefficient of the probe laser 
as the laser frequency is tuned. The input for the sim- 
ulation is the laser detuning, which partially determines 
the steady state of the electron-hole system. The output 
is the expectation value of the absorption coefficient for 
the probe laser. Since the absorption of the probe laser is 
proportional to the imaginary part of the density matrix 



= lm\(x — 
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as the absorption coefficient of the probe laser. The key 
point is to calculate the expectation value of the nu- 
clear spin Overhauser field, since the actual detunings 
for the optical transitions take into account the nuclear 
spin Overhauser field, which shifts the electron spin Zee- 
man energy. As for the input parameters, 7 S is in prin- 
ciple determined by the fluctuation of h. However, as 
mentioned in Sec. IIV Al we neglect the finite width of 
Ps(si„,s As ) when studying the second type of observa- 
tions. This means that we have lost the information on 
the evolution of the nuclear spin fluctuation. As a re- 
sult, we apply a constant electron spin decoherence rate 
7 S taken from the estimation in Ref. [l3| for each set of 
chosen parameters. 



1. Narrowed distribution of the nuclear spin field 

Following Ref. Hi], we assume that the system contains 
only indium nuclear spins, since the momentum of an in- 
dium nuclear spin is three times that of an arsenic nuclear 
spin, and A\ n > Aas- In this case, we rescale the indium 
nuclear spin hyperfine constant as shown in Appendix [Cl 
Equation (fTUj) becomes 
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As in Refs. [3l|, [53J, and [56J, the solution for the steady 
distribution of nuclear spin polarization is 



p( ss )(s la )=Afexp 



A„K) 



ds' 
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where J\f is the normalization factor. From pi ss \s Ia ) 
we can evaluate the steady distribution of the nuclear 
spin Overhauser field p^ s \h). Using the same param- 
eters (except Ai :0 ) in Fig. [3J we plot p^ (h) for dif- 
ferent nominal detunings Ai,o = —0.3, 0, 0.675 and 
±2 GHz in Fig. [3fa), (b), (c) and (d), respectively. For 
comparison, the dashed curves denote the distribution 
of the nuclear spin Overhauser field in thermal equilib- 
rium (the mean of the Overhauser field is set to be zero, 
shown in Appendix [G]) . One can see that for each case, 
the steady distribution of the nuclear spin Overhauser 
field is narrowed compared to that in thermal equilib- 
rium. Also, the largest nuclear spin polarization in these 
four cases occur in Fig. [3Jd) , where the mean of s ln is 
- T0.067 for Ai, = ±2 GHz, fulfilling the condition 
|s In | 1 that is needed for deriving the Fokkcr-Planck 

equation. An n-timcs narrower p^(h) indicates n-times 
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FIG. 4. (Color online) (a) Nuclear spin Overhauser field, 
and (b) Xi as a function of detuning Ai^o when we change 
the laser detuning by 0.04 GHz at 4 s interval. Data shown 
with solid (dash-dot) curves correspond to changing Ai,o from 
left (right) to right (left). Here we use 7d cp = 0.2/s. Other 
parameters except Ai,o are the same as in Fig. [3] 



smaller an electron spin decoherence rate^ Specifically, 
Ai.o = 0.675 GHz in Fig. |3tc) is equal to 172/2, cor- 
responding to the detuning of the second pump in the 
two-pump experiment shown in Fig. 3c of Ref. [l3|, where 
a significant enhancement of the electron spin coherence 
is observed. Besides a narrowed distribution of the nu- 
clear spin Overhauser field, we also note that there are 
two separated peaks in p£ B \h) when Ai.o = 0, indicating 
bistability of the nuclear spin polarizationi 4 i 8 i 57 Finally, 
we note that p^ s \h) for lA^ol and — |Ai o| are symmet- 
rical to each other about the line h = 0, as shown in 
Fig. EUd) for illustration. 



2. Laser spectroscopy experiments 

Here, we study systems with two different sets of pa- 
rameters and show that, as they were captured, three 
phenomena observed experimentally in Ref. [TU The first 
case is for showing that the absorption curve will shift 
when the probe laser frequency is changed step by step. 
In the second case, we will show that (i) the line-shapes 
of Xi f° r forward and backward scanning are not symmet- 
ric to each other when the pump laser is off-resonance, 
and (ii) there can be an abrupt switching of the probe 
laser absorption when we stop changing the laser fre- 
quency at a rising edge of the absorption curve. We note 
that these three phenomena agree with the respective 
experiments (see Fig. 1, Fig. 2a, and Fig. 2c of Ref. Ii~3l . 



FIG. 5. (Color online) Same as Fig. [4] Parameters are 
the same as in Fig. [3] except 7 S = 0.22, A2,o = —0.8, Q.2 = 
0.7 (unit: GHz), and 7 dep = 0.02/s. 



respectively). 

For the first case, the parameters A 2i0 , I\ 7 S! fii, fl 2 
are set as 0, 0.4, 0.06, 0.24, 1.35 GHz, respectively. We 
set 7dcp = 0.2/s here. We change the central frequency 
u>i of the probe laser by 40 MHz at 4 s interval, from high 
to low and from low to high, corresponding to chang- 
ing Ai t o from negative to positive and from positive to 
negative, respectively. By using Eq. (fTTj) . the calculated 
Overhauser field h and the absorption coefficient \i °f 
the probe laser is shown in Fig.fJta) and (b), respectively. 
The data shown with a solid (dash-dot) curve correspond 
to changing A^o from left to right (from right to left), 
respectively. One can see that the absorption curve shifts 
towards the direction of changing Ai^- This result agrees 
with the experimental observation in Ref. K3l 

For the second case, we choose a detuned pump laser 
with A2 ) o = —0.8 GHz (this does not mean the detun- 
ing in Fig. 2a of Ref. HH is negative. We use detuning as 
the x-axis, while Fig. 2a of Ref. fl3l uses laser frequency), 
and choose the Rabi frequency of the pump laser O2 as 
0.7 GHz. Since the pump laser is weaker than the one 
used in Fig. |4j we set 7 S = 0.22 GHz, as estimated from 
Ref. [HI (see Fig. 1 of its supplementary information). 
Also, we set 7d ep = 0.02/s, which is much smaller than 
the depolarization rate used in Fig. |4j Interestingly, we 
note that in the supplementary material of Ref. Il3l . a 
smaller (comparing with the one used for the resonant 
pump laser) nuclear spin depolarization rate for a de- 
tuned pump laser is also used for numerical simulation. 
We have chosen the parameters different from those in 
Fig. JH in order to reproduce the experimental results in 
Ref.ll3l. Indeed, when we calculate the nuclear spin Over- 
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FIG. 6. (Color online) Left and right panels of (a): h and Xi as 
a function of t after stopping changing Ai,o from 3.5 GHz step 
by step [The step is 0.04 GHz as in Fig.E] to Ai, = 2.3 GHz, 
denoted by the arrow in Fig.(5jb). (b), (c) and (d) are similar 
to (a) except that we stop changing A^o at 2.335, 2.34 and 
2.35 GHz, respectively. The inset of the right panel of (d) 
shows the longer time behavior of \i- One can see that the 
absorption coefficient has a switching behavior in (b) and (c). 



hauser field and the absorption coefficient of the probe 
laser shown in Fig. [5J we can see that the line-shapes of 
Xi for forward and backward scanning are not symmetric 
to each other. 

In the second case, we also study a switching behavior 
of the probe laser absorption. As indicated by the arrow 
in Fig.[5]Jb), we now change A^o from right to left step by 
step until reaching A^o = 2.3 GHz, where we still keep 
both lasers on. At the instant we stop changing A 10 , we 
begin to record the laser absorption as a function of time, 
we present h and Xi as a function of time in the left and 
right panel of Fig. [H^a), respectively, where we find no 
abrupt switching of Xi ■ We then repeat the same calcu- 
lation as in Fig. Uta) when the point of Ai.o at which we 
stop changing Ai.o is set as 2.335, 2.34 or 2.35 GHz. The 
respective results are shown in Fig.H^b), (c) and (d). We 
can see that there is an obvious switching behavior of Xi 
at around t = 310 s (or t = 610 s) after we stop scan- 
ning the probe laser in Fig. EJb) [or (c)]. The switching 
behavior here qualitatively agrees with the experimental 
observation in Fig. 2c Ref. Il3l . Interestingly, our simula- 
tion reveals that there would be no such behavior when 
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FIG. 7. (Color online) vi a /Ti ntP as a function of Si n . Ai.o = 
—0.3,0,0.675 and 2 GHz in (a), (b), (c) and (d), respectively. 
Other parameters are the same as in Fig. [3] In each subfig- 
ure, the right panel zooms in the left panel, showing the cross- 
ings [which gives the position of a local maximum of pi ss ^ (s In )] 
between the curve ui n /ri n>p and the horizontal dashed line. 
One can see that there are two such crossings in each panel, 
denoted by the circles. 



we stop shifting the probe laser frequency too far from 
the rising edge, say, at Ai = 2.35 GHz, as shown in 
Fig.Efd). 

The hysteresis in Figs. @] and [5] and the switching effects 
in Fig. [6] result from the nonlinear feedback between the 
electron-hole system and the nuclear spin bath, where the 
feedback is controlled by the hole spin-nuclear spin non- 
collincar hyperfine interaction. In Appendix A, we show 
that the spin flip-flops due to the hyperfine interaction 
may be neglected for strong magnetic field. However, 
the feedback is much less effective in polarizing the nu- 
clear spins if the external magnetic field is lower than, 
for instance, 0.1 T. This is because the function s In(Aa):0 
in Eq. (|B9j) that is responsible for the nonlinear feedback 
will vanish as ~ thus causing the hysteresis and 

switching effects through the noncollinear interaction in 
our model to vanish. In such a case, the dynamics of the 
nuclear spins will be controlled primarily by the hyperfine 
flip-flops. It has also been noted that similar hystere- 
sis of nuclear spin polarization in an optically pumped 
QD could be explained by exploiting the spin flip-flops 
through the electron spin-nuclear spin hyperfine interac- 
tion [see, e.g., Refs.lilMflil and [53] ■ 
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FIG. 8. (Color online) Same as Fig. [7] except that here f2i 
1.35 GHz and fi 2 = 0.24 GHz. 



FIG. 10. (Color online) Same as Fig. [7] except that here 
fii = fl 2 = 1.35 GHz. 




FIG. 9. (Color online) Same as Fig. [7] except that here f2i 
fi 2 = 0.24 GHz. 



V. LARGE NUCLEAR SPIN POLARIZATION 

Now, wc turn to the cases in which the nuclear spin 
dynamics controlled by Eq. ([3]) can induce significant 
nuclear spin polarization, where we cannot apply the 
Fokkcr-Planck equation from the preceding section since 
its derivation requires |5 In ( As j| 1. In order to simplify 
the discussion, wc follow Sec. lIVB H and assume a system 



with only indium nuclear spins. 



A. When |a In | < 1 is violated? 

Before introducing a condition where |s In | 1 is not 
fulfilled, we first study cases where |s In | <C 1 is satisfied 
and Eq. (fT4|) determines the evolution of p s {s Ia ). For 
this case, it is pointed ouiM that a peak of pi ss (sm) can 
usually be approximated by a Gaussian. The position 
Sin = Si n>m of such a Gaussian, i.e., a local maximum of 
=S '(sin) in Eq. (|16|) is given by 



P 



«ta(Sln,m) _ q 9 V ln (s[J 



D ln (s ln , m ) ' S<A„(s' In ) 
From Eq. ^ we have D In (s ln ) ~ (J+1)ri " 



< 0. 



(17) 



3NI 



hence we 



can study Ui n /ri n)P to find the solution of Eq. (fl7|) . Wc 
plot ?;in/ri ni p for four different nominal dctunings Ai ; o = 
-0.3,0,0.675 and 2 GHz in Fig. [7[a)-(d), i.e., the same 
as those in Fig. EJa)-(d). From Fig. [3 one finds that 
there are two solutions (denoted by the small circles) to 
Eq. (fTTj) for each case, hence there should be two peaks 

of pi BS \s ln ). Nevertheless, we note that for the two local 

maxima of pi ss \s ln ) for Ai o = 0.675 or 2 GHz, one peak 
is much lower than the other peak in Fig. [3fc) or (d) that 
it is not visible in Fig. [3] [mathematically, taking Fig.[3Jd) 
as an example, for the two deformed 'triangles' formed 
by the solid curve and the horizontal dashed line, the 
area of the triangle pointing down is larger than that of 

the one which points up]. All local maxima of pi ss ^ (s In ) 
given by Fig.[3]occur with |s In | <C 1, indicating that only 
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small nuclear spin polarization can be generated. This 
is consistent with the condition for deriving the Fokker- 
Planck equation. 

Now we consider the case where |s In | 1 is not ful- 
filled. As an example, we exchange the intensities of the 
pump and probe lasers. If we assume \s ln \ <C 1, we can 
still derive the same Fokker-Planck equation Eq. (fTDj) . In 
this formalism, we can also assume that the position for 
a local maximum of Ps (s In ) is given by a solution to 
Eq. (JT7J). Again, we plot w In (s Ia )/r Ini p in Fig. [5] One 
finds that there is always a solution to Eq. (TIT)) at s In ~ 1. 
If one naively tries to calculate pi^ (s In ) in any of the four 
cases in Fig. [8j the peak ~ 1 dominates and the mean of 
s In is almost 1. On the other hand, if we use Eq. (TIT]) 
to study the evolution of the mean of nuclear spin polar- 
ization, strong polarization results. This contradicts the 
condition |s In | <C 1. 

When we set fii = fl 2 and choose weak pumping, e.g., 
fii = Q 2 = 0.24 GHz, the solution to Eq. (JT7|) at s tn ~ 
—0.2, as seen from Fig.^d), is already much farther from 

than the corresponding solution shown in Fig. UJd). 
When both lasers are strong, e.g., f?i = CI2 = 1-35 GHz, 
Fig. rrU] shows that the solutions to Eq. (fTTj) can be even 
further from zero in Fig. HOf c) and (d), thus indicating 
the failure of the Fokker-Planck equation. 

In conclusion, the dynamic nuclear spin polarization in 
a system illustrated in Fig. [1] will obey (violate) |s In(Aa) | <C 

1 when Qi -c ^2 (^1 S> fl 2 ) and A 2 ,o = 0. In any 
other case the condition |s In(As) | <C 1 may fail to hold, 
depending on the Rabi frequencies and the dctunings of 
the two laser beams. 



B. Numerical study when |si„ ~ 1 

Here, the condition for the derivation of a Fokker- 
Planck equation breaks down, but we can still solve 
Eq. ([3]) for the dynamics of the nuclear spin Overhauser 
field when \s In \ ~ 1. We study an artificial QD with 
only 40 nuclear spin J = | by the Green's function 
Monte Carlo (GFMC) .21^ We assume that there are 
20 gallium and 20 arsenic nuclear spins. The gyromag- 
netic ratio gjfi K — 0.0093 GHz/T and hypcrfinc constant 



10.6 GHz, are estimated by taking the average 



among 69 Ga, 71 Ga, 75 As (the natural abundances for 
these three are 60.1%, 39.9% and 1, respectively^) for 
the gyromagnetic ratios and hyperfine constants, where 
the hyperfine constants for gallium and arsenic nuclear 
spins are taken from Refs. |3 andl52l 

A detailed review of the numerical method used here 
can be found in Ref. [59|. In order to implement GFMC 
simulation, we convert Eq. ((3]) into the form of a first- 
order differential equation, 



dt 



p n {t) = Ap n (t), 



(18) 



d = (2J + l) 40 different states. The d x d matrix A 
determines the dynamics of the nuclear spin population. 
From Eq. ([3]), one can prove that each off-diagonal el- 
ement of A is non-negative. Nonetheless, the diagonal 
element Akk is given by 



.4 



kk 



5> 2 + j if - j?)w. 



./.+ 



^(J 2 + J - Jf + Jf)Wj,- 



k), (19) 



where p n (t) is a vector, each of whose elements repre- 
senting the population of nuclear spins on one of the 



which is non-positive. Here, Wj t ± is W aj ,± defined in 
Eq. dU, and \k) represents one of the d nuclear spin 
states. Most of the eigenvalues of the matrix A are nega- 
tive (except one being equal to zero, which we tested by 
exactly diagonizing A for smaller systems). The eigen- 
vector of A corresponding to the biggest eigenvalue of 
A gives the steady population of the nuclear spins. In 
GFMC, we can study properties of the eigenvector cor- 
responding to the biggest eigenvalue of the matrix A. 
To do this, we first choose a positive constant and add 
this constant to every diagonal element of A, shifting its 
spectrum up by making each of its diagonal matrix el- 
ements positive. This docs not alter the result of the 
GFMC simulation^ At the start of the simulation, we 
randomly set the state of each nuclear spin, which is 
equivalent to assuming that the temperature for the nu- 
clear spins is infinitely high. During the simulation, we 
run M independent sets of simulations simultaneously for 
im Monte Carlo sampling. After each step of sampling, 
we record the distribution p s (s) of the nuclear spin polar- 
ization s = -jfj J2j(Jj) calculated from the M different 
nuclear spin states. After a certain step i' M the nuclear 
spin ensemble of these M states reaches equilibrium, and 
the distribution of s is calculated by analyzing the data 
from the i' M step to the im step. 

In order to test the GFMC in our model, we can first 
consider the case f?i(2) = 0.24(1.35) GHz, in which case 
the nuclear spins polarize weakly when N = 9500, as 
shown in Fig. [JJ We assume that the condition \s\ -C 1 
holds when we consider a small system with only N = 40 
nuclear spins. In this case, we can use both the Fokker- 
Planck equation and the Monte Carlo simulation to cal- 
culate the steady distribution of the nuclear spin polar- 
ization s. The reason that we may still use the Fokker- 
Planck equation is that the condition A-j <C 1 for de- 
riving the Fokker-Planck equation (see Appendix [Bj still 
holds when N = 40. In Fig.[TTJ we present the steady dis- 
tribution pi s (s) of s for a system containing 40 spin-| 
for Ai t o = —1 and 1 GHz in (a) and (b), respectively. We 
calculate Ps (s) from both GFMC and Fokker-Planck 
equation. For the former, we set M = 200, Im = 10 5 and 
i' M = 10 4 in the GFMC study. We also use Eq. ^ to 

calculate pi ss ^ (s) for a comparison. From Fig. HU one 
can see that the results from GFMC and Fokker-Planck 
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FIG. 11. (Color online) Steady distribution of the nuclear 
spin polarization s for a system with 40 spin-|, calculated 
from the Green's function Monte Carlo simulation (solid curve 
with circles) and Fokker-Planck equation (dashed curve), re- 
spectively. Here the nuclear spin gyromagnetic ratio is <7j/x N = 
0.0093 GHz/T, and the hyperfme constant is A 3 = 10.6 GHz. 
A^o = —1 and 1 GHz in (a) and (b), respectively. Other pa- 
rameters are the same as in Fig. [2] For comparison, we have 
shown the result for a thermal nuclear spin state by dash-dot 
curve. 



equation agree very well with each other. Another fea- 
ture in Fig.[TT]is that the distributions of s for A 10 = — 1 
and 1 GHz are symmetric to each about the line s = 0. 

We then consider the case Oip) = 1-35 (0.24) GHz. 
From Fig. [5] [the curve v Ia /Ti n!P for the present case is 
different from that in Fig. [51 yet also indicates that s«l 
does not hold], we know that the criterion for using the 
Fokker-Planck equation fails to hold. But the Monte 
Carlo simulation is still valid. We sample M = 200 
simulations simultaneously. However, we found that the 
state represented by these samples does not reach equi- 
librium [i.e., p s (s) is going up or down for a given s if 
we continue the sampling] when we run only i' M = 10 4 
simulations, as above. So we set %m = 10 6 here. Now s is 
analyzed by the data generated from the i' M = 91 x 10 4 th 
to the ijvfth step. In this case, we found that the Monte 
Carlo simulation reaches equilibrium (we do not mean 
that at least 91 x 10 4 sampling is needed to reach equilib- 
rium). The distribution of s from this Monte Carlo simu- 
lation is shown in Fig. Q21 One can see that the distribu- 
tion of the nuclear spin polarization is peaked at s « ±0.5 
when A^o = ±1 GHz, and it is narrower compared to 
the thermal nuclear spin state without optical pumping. 
We note that the distributions of s for Ai.o = — 1 and 
1 GHz arc almost symmetric (except each of them has a 
small kink on the left shoulder) to each other about the 
line s = 0, similar to that in Fig.[TT] Because the peak of 
the nuclear spin polarization distribution is shifted from 
to ~ ±0.5, we conclude that strong nuclear spin po- 
larization can be generated by optical pumping on single 
QDs. 
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FIG. 12. (Color online) The solid curves with circles and the 
dash-dot curves mean the same as the counterparts in Fig. [TT] 
except that they are calculated with fii — 1.35 GHz and 
Q 2 = 0.24 GHz. 



VI. CONCLUSIONS 

We have studied dynamic nuclear spin polarization in 
a strained self-assembled QD pumped by two narrow- 
linewith continuous wave lasers, where the noncollinear 
hole spin-nuclear spin hypcrhnc interaction is assumed 
to be responsible for nuclear spin polarization. We found 
that the nuclear spins can be polarized to a degree that 
is either small or large, depending on the intensities and 
central frequencies of the lasers. 

We first study the experiments in Ref. [TH, where only 
weak nuclear spin polarizations were observed. In this 
case, we derive a Fokker-Planck equation for the time 
evolution of the probability of nuclear spin polarization, 
based on which wc show that the distribution of the nu- 
clear spin Ovcrhauser field can be narrowed, as well as 
other phenomena, including that the absorption curve of 
the probe laser shifts when we change the probe laser 
frequency. When large nuclear spin polarization can be 
generated by optical pumping, the condition to derive a 
Fokker-Planck equation breaks down. Then we use the 
GFMC simulation to directly solve the steady distribu- 
tion of the nuclear spin Ovcrhauser field. Indeed, we find 
that a large nuclear spin polarization up to 50% can be 
generated for a small 'QD' containing 40 spin-|. Note 



10 



that a 50% nuclear spin polarization is close to the large 
nuclear spin polarization observed in experiments by op- 
tical pumping 2 - ' 11 ! 20 ' 36 or other controls^ in QDs. 

Throughout this paper we have assumed that the dy- 
namic nuclear spin polarization results from the hole 
spin-nuclear spin noncollinear hyperfine interaction in a 
self-assembled QD under optical pumping. When the 
quadrupole interaction^ between the electron and nu- 
clear spins is strong compared to the Fermi-contact hy- 
pcrhnc interaction in the QD, there will be an effective 
noncollinear interactio n 16 ! 26 between the electron spin 
and nuclear spins which can play a similar role in po- 
larizing the nuclear spins as studied in this paper. We 
note that our method is applicable when nuclear spin 
polarization is mainly controlled by the effective electron 
spin-nuclear spin noncollinear interaction. 
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Appendix A: Steady state of the electron-hole 
system 

This appendix gives the steady-state solution to the 
reduced density matrix of the electron-hole system in an 
adiabatic nuclear spin Overhauser field h. This steady 
state is determined by H c h and the longitudinal part of 
the hyperfine interaction H nl in Eq. (p}. 

First of all, it is helpful to use a rotating frame to show 
that the optically pumped electron-hole system in Fig. Q] 
is an effective three-level system. The Hamiltonian for 
the electron-hole system includes three parts, namely the 
Zeeman energy -f/eh-b m an external magnetic field, the 
control H c h-\ by the electromagnetic fields and the chan- 
nel of photon exchange -ff G h- q with the electromagnetic 
vacuum, 

Heh = -ffoh-b + i?oh-l + -ffeh-q- (Al) 



H ch . h - E x \x+)(x + \ + {E x + g e ^B)\x-){x-\+E T \T+){T+\ + {E T + g hf x B B)\T-){T - \, (A2) 

where E x and Et are the eigenenergies of the eigenstates |a;+) and \T+). For the system in Fig.[TJ the semiclassical 
Hamiltonian of the control by the electromagnetic fields on the electron-hole system in dipole approximation^ is 

tfch-i = ^(e 1 " 1 ' + e-^x-) + ^(e*""* + e -^ t )\x+) S j <T - | 

+ (^(e Mlt + e~ i ^ t )\x+) + ^( e ^ 4 + e^ 2 *)^-)) (T+\+ H.c, (A3) 

where u>i(2) is the central frequency of the two coherent laser beams. Here, we assume that the matrix element of the 
electric dipole momenlj 6 - 2 - e(x + \x(y)\T±) is equal to e(x — \x(y)\T^f) (e is the elementary charge), hence the Rabi 
frequency for the transition |x=f) — > \T+) is equal to that for |x±) — > \T— ). Because the two laser beams are almost 
resonant with the two transitions \x±) —> \T—), we use a rotating frame with 

R = (E T + g h ^ B B - U!)\x-){x - \ + (E T + g h ^ B B - uj 2 )\x+)(x + \ + (E T + g hl i B B)\T-){T - \ 

+E T \T+){T+\, (A4) 

to eliminate the obvious time dependence in H e h for the transition |x±) — > \T— )(we use the subscript r to denote the 
rotating frame), 

« -Ai >0 |x->(x - | - A 2 ,o\x+)(x + I + (^->( T - I + ^-\*+)(T - | + H.c) + e mt H c ^ E - lR \ (A5) 
Ai.o = E T + g h fi B B - (E x + g c p B B + wi), A 2 , = E T + g h fi B B - (E x + w 2 ), 



where we have ignored the rapidly oscillating terms with 
phase term exp[±i(wi + 0J2)t] or exp{±i(g c ± gh)n B Bt] 
since Wi( 2 ) ~ 3 x 10 5 GHz^i 3 - and \(g c ± g\ l )^ B B\ > 



13 GHz [with B = -2.64 T and g c(h) = 0.49 (-0.13)], 
much larger than | A 1( - 2 ) | used in the main text and other 
parameters in the Hamiltonian. From Eq. (|A5[) . we can 
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see that the effective electron-hole system involves \x±) 
and \T— ), hence can be called a A system (AS). 

When we apply the rotating frame transformation 
to the transverse part of H m , we have acjSflJ — > 

a CJ 5 c ± ^e ± '" (A2 ' " Ali0 " Sc, ' BB) - Here g c p, B B w 18 GHz, 
hence e ±it(A.2,o-Ai, -9e/iBS) is rapidly oscillating when 
Ai.o and A 2i o are set as in the main text. So the 
transverse part of the electron spin-nuclear spin hy- 
perfine interaction can be neglected. We note that 
a Schrieffer- Wolff transformatio n 63 ! 64 on the trans- 
verse part of H H1 may give rise to a nonlinear term 
"e S j j'^tj 2g tmB J' ' wmcn a ^ s0 contributes to nu- 
clear spin dynamics 4£ However, one can evaluate 6 ^ and 
show that the pumping of nuclear spins through this non- 
linear term is much slower than that through the non- 
collinear term ~ S^lj for the specific system studied in 
this paper. For such case, we shall assume that the trans- 
verse part of electron spin-nuclear spin hyperfine interac- 
tion can be neglected. Finally, since ah.j/a CJ ~ —0.1 (see 
Refs. Hi] and l5ll ) and the strength of the coupling be- 
tween |T+) and \T— ) given by a mean field treatment of 
2 I^K,j §yfx sm< 5 i s wea k compared to the Rabi frequen- 
cies and decay rate of the trion in this paper, we further 
neglect the transverse part of the hole spin-nuclear spin 
hyperfine interaction. Equation (J2J then becomes 



Hr, 



E 

3 



QX TX 

>,3 D e 1 j 



\ a h,j 



V3 



S?(77 cos 5 



I] sin 6) 



(A6) 



where we have used 1 + |/3| 2 sa 1 when \/3\ = 0.2. 

We perform the mean field approximation on the lon- 
gitudinal part of the electron (hole) spin-nuclear spin hy- 
perfine interaction. This gives an effective nuclear spin 
field, i.e., the Overhauser field that shifts the Zeeman 
energy of the electron (hole) spin, 



3 



(A7) 
(A8) 



i/a e 



-0.1, 



Since ~ 0.2 (taken from Ref.[l3l) and Oh, 
we neglect but keep the nuclear spin Overhauser field 
h via the electron spin-nuclear spin hyperfine interaction. 
As a result, the actual detuning between the photon en- 

I 



ergy of the laser and the energy cost in relevant transi- 
tions is 



Ai = Ai n 



h 

2' 



^2,0 



h 
2' 



(A9) 



Next, we solve the steady-state reduced density matrix 
of the AS in an adiabatic nuclear spin Overhauser field h. 
The last term in Eq. (|A5[) induces spontaneous decay of 
the trion level, which we incorporate in the optical Bloch 
equation in the Lindblad form^ 6 - 

^ = -i[HAS,r,PA] + \GkPA& k - \ {6{G k ,p A } 

(A10) 



k=l 



+7«£ G k pAGl-±{G{G k ,p A } 

k=3 L 



wherc {a, b} = ab + ba, and 
# AS>r = -Ai|a?-)(a5 - I - A 2 |a;+)(.T 



^\x-)(T-\ + ^-\ X +)(T-\+R.c 



Gi, 2 = \x±)(T - |, G 3 ,4= \x±){x±\, 



(All) 



and the A-system reduced density matrix p A in the basis 
constructed from the vectors \T— ), \x+) and \x—) is 



Pt.t Pt,x Pt.x 

PA = I Px,T Px,x Px,x 
PxT Px,x Px,x 



(A12) 



The condition Tip\ = 1 allows us to eliminate one matrix 
element, say, ps,x\x— ){x — | in p\ and rearrange its other 
eight matrix elements into one column, 

PA = (Pt,T' Px,x, Px^Ti Px,T> Px,x, Pr,xi Pt,x.i Px.x) 1 



where T denotes transposing a matrix. Then Eq. (IA10[) 
becomes, 



*^a(*) =Mp A {t) + X. 



Here, 



and 



x = (o, 0, 0, 



T' °' °' T' 0) ' 



(A13) 



(A14) 
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(A15) 



Ai - ij s } 
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where Ti = 2T is the relaxation rate of the trion state, 
r 2 = r + is the decay rate of the coherence between 
the trion and one of the electron spin eigenstates, T is 
the spontaneous decay rate from the trion to each of the 
two electron spin states, and 7 S is the energy-conserving 
dephasing rate of the electron spin in the presence of fluc- 
tuating nuclear spins. Note that we have ignored the elec- 
tron spin relaxation process which is much slower than 
the electron spin decoherence^i 

In the steady state of the AS, we obtain 

Mp A (t -> +oc) + X = 0, (A16) 
I 



whose solution, together with the condition Trp\ = 1, 
gives us the reduced density matrix for the steady 
state of the AS. 

Appendix B: Fokker-Planck equation 

In this appendix, we derive the Fokker-Planck equation 
starting from Eq. (J3J) . For the derivation of an analogous 
equation involving only one species of nuclear spins, see 
Ref. IsH We first multiply the Kronecker delta function 
^Sin,si„^sA 3 ,sA B on both sides of Eq. ([3]), and then trace 
over the nuclear spin degrees of freedom in it, giving 



EM 



IJW t 



P{t) 



^Si n ,Si n <5s As ,s Aa 



if, ijw aj ,-P{t) 



^Sln-Sln^SAs, 



■}■ 



(Bl) 



where Tr denotes tracing over the nuclear spin degrees of freedom. Using Eq. the left hand side of Eq. (|B1[) 
becomes 4zp s (s lxl , s Ab ) [we suppress the variable t in p s for brevity]. Taking the sum in the first commutator on the 



right hand side of Eq. (|B1|) as an example, the trace evaluates to 



-E Tr {[^' lfW aj , + P(tj\ 'K,. , SA } = -EE( fc l{[V- ti~W aj ,+P(tj\ 6 Sltt , sl J SAe , SAs } \k), (B2) 

3 3 k 

where k labels a specific nuclear spin state and runs over the whole Hilbert space of the nuclear spin states. 



First, we evaluate the sum for Ij IjW ajt +P(t) on the right hand side of Eq. (|B2|) : 

- EE( fc i { [^^w aj , + p(t)] s Sln , sl j SAe , SAB } \k) = EEw l 1 ^ + !) - ! f - w at ,+p(t)s h 



3 k 



3 k 

xSs Ae ,s AB \k). 



(B3) 



Below, we label Ii n = I = | and /as = J — § unless otherwise specified. Now assume that the following approximation 
is valid when both polarizations of the indium and arsenic nuclear spins are small^i 



E (^) 2 

\j'Gln(As) 



NI lnlAs) (I la(As) + 1) 



(B4) 



where (...) denotes the respective expectation value, then Eq. (|B3[) becomes, 

~2NI a (I a + l) 



3 k 



EEw { [ttw^+p®] 5 haiS1 j §As , SA \\k) « E 

- ,As 

As 5 S As I 

~2NI a (I a + l) 



k a— In, As 

X ^5ln,Sln^A S ,SA S 1^) 



E 

a— In, As 



W a ,+P(t) 



Wa,+ (Si n , S As )p s (s ln , S A s), 

(B5) 



where s a , with a = In (As), is defined in Eq. (jB]). 
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Next, we evaluate the remaining sum on the right hand side of Eq. ()B2|) : 



3 k 



^si n ,si I1 ^SAs,SAi 



4- (similar term with sum over As) 
2(7 + 1) 



,Sl n — O^Sab,SAs 

I*) 



A/7 



3 

2(J + 1) 



Sin + a 



Wi„ i+ (s In - a, s AB )p s (s In - a, s Aa ) 



s As + b 



W / As,+ (s In , s As - b)p s (s ln , s As - b), 



(B6) 



where a = -^y , and & = -^j . 

Similar to Eqs. (|B5[) and (|B6|) . the sum over the second commutator on the right hand side of Eq. (|B1|) gives, 



^il D ,«In^ As ,SAs} = - E NIa 
a— In, As 

-2(1+1) 



2(I a + 1) 



+NI 
+NJ 



2(J + 1) 



s As + 6 



fin; 

Win,- (s : „ + a, s Ae )p s (s ln + a, s As ) 
WA B ,-(si n5 s As + b)p s (s la , s As + b) 



(B7) 



To simplify Eq. (|B7[) . we perform series expansion up to the second order in a and b, and obtain 



+jVJ&W As -(sid SA B )ft(si„, s As ). 



+NJ 



2(1 + 1) 



2(J + 1) 



s As + 6 



ds In 2 9s 2 



2 ds 2 . 



Defining 



rin(As),p(Sin, S As ) = Wl n (As), + (Sin ; S As ) + Win(As),- ( s in, S A8 ), 
Si„(As),o(Si„, S As ) = [Win(As), + (Stn, S As ) - Wln(As),-(Sin,SAa)J /rin(As),p(Sin, S A s), 

and using Eqs. ()B5|) . (|B8|) and series expansion from Eq. (|B6[) [similar to Eq. (|B8[) ] . Eq. (|B1[) becomes, 



(B8) 



(B9) 



E K 



a— In, As 



,pPs 



2(I a + 1) 



1 d 2 T aiP p s ds a ^ a , P p s 



2NI a ds 2 a 



8S C 



1 d 2 s afi r atP p s dT atP p s 



2NI a dsl 



ds a 



E 



=In,As 



1 ds afi r a ^ p p s 

NI a ds a 



E 



=In,As 



d_ 

8s a 



_d_ 

ds a 



D a p s (s 

In: In) S Ab 



(BIO) 



where we have dropped terms ~ 1/N 2 . The drift and diffusion coefficients in Eq. (|B10[) are given bj*2i — 



fln(Afi) 



D 



In(As) 



In(As),p 



In(As),p 



2(/ln(As) + 1) 

Sln(As) 5 Sln(As),0 



2NL 



In(As) 



2(7ln(As) + 1) 



S In(As)Sl„(As),0 



(Bll) 
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In the definition of v In{Aa} and -Dwas) m Bq. (|B11[) . the factor 2 ( /ln (* s ) +1 ) comes from the sum of individual nuclear 
spin fluctuations in Eq. (|B4|) . and r In (A s ), P is the rate of nuclear spin flip due to the noncollinear interaction in our 
model. If ri n (As),p is large, the distribution of the nuclear spin polarization evolves rapidly. The nonlinear function 

•Sma 3 ),o plays a key role in our feedback loop since the solution to s In(As) — 2 ( f M^") +1 ) s Iii(Ab)|0 = gives us the stable 
nuclear spin polarization (see Sec. IV A| . Because |s In ( AB) s In(Aa); o| ^ F the diffusion coefficient can be approximated as 

n _ r in(As),p 2(J In(As) +1) 

Mn(As) ~ TTTTT 5 - OlnfAsVth 1 In(As),p, \P^) 



where oi n (As),th = y i^Nh]^ 1S ^ e standard deviation of nuclear spin polarization distribution for N indium (arsenic) 
nuclear spins at infinite temperature ! 31 ! 52 The appearance of 0W.As),th i n the definition of Di n (As) is not surprising 
since -Din(As) is associated with the diffusion process.— 

Defining the means of the nuclear spin polarization s ln and s Aa by 

s ln p a ds In ds Ae , 

s As p s ds Ia ds As , (B13) 

and further assuming that the probability for the nuclear spins to be totally polarized is negligible, i.e., p s (si n 
• I - s As ) = p s (s In , s As = ±1) = 0, we obtain 3 ^ 

| 5 ,. = //», J ».( 8 ,„ 8 „)*,^, 

^s As = I f v As p s (s Aa ,s Aa )ds ln ds Aa . (B14) 



Appendix C: Inhomogeneous broadening 



The nuclear spin Overhauser field obeys a Gaus- 
sian distribution when the nuclear spins are in thermal 
equilibrium. ; 52 ! 68 



Ph(h) 



-(h-hf/(2r* 2 ) 



2tiT* 



(CI) 



h = Ainlsn, + Aa s Js As , 

,2 ^ A 2 ln I(I + l) + A 2 As J(J + l) 

2 3A 



when |s r n(As)| "C 1 for temperature ~ 5 K in typical 
experiments! 13 ' 14 Here h and are the averages of h 
and s a , respectively. The electron spin decoherence L(t) 
is dominated by the inhomogeneous broadening of the 



nuclear spin Overhauser field ^ 

L(t) = ( p h (h)e~' l ^^ BB+h ^dh. 



(C2) 



pure Gaussian decay^ 



For a thermal nuclear spin bath of Eq. (|C1[) , we have a 

^-i(g c ^ B B+h)t-(t/T;) 2 



L(t) 



(C3) 



where T 2 * = V^/rj- Both indium and arsenic nuclear 
spins contribute to in a QD with N indium and ./V 
arsenic nuclear spins. Under the condition of |s In(As) | <C 1 
and further ignoring the arsenic nuclear spins in the QD, 
it is useful to rescale the hyperfine constant of indium 
nuclear spins so that the new system with only N indium 
nuclear spins gives a F^ that is the same as in the original 
InAs system. This is achieved by writing 

A' 2 In I(I + 1) = A 2 In I(I + 1) + A\ a J{J + 1), (C4) 

which gives a 5% increase in the hyperfine constant for 
indium nuclear spins. However, this does not mean that 
the new system with only N indium nuclear spins retains 
all the properties of the initial system. 
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